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^! Abstract 

o 

The coordinate transformation method is apphed to bounded domains to design metamaterial 



devices for steering spatiaUy confined electromagnetic fields. Both waveguide and free-space beam 

applications are considered as these are analogous within the present approach. In particular, we 

^~»| describe devices that bend the propagation direction and squeeze confined electromagnetic fields. 



Two approaches in non-magnetic realization of these structures are examined. The first is based on 



CN ' using a reduced set of material parameters, and the second on finding non- magnetic transformation 

>, 

^^ , media. It is shown that transverse-magnetic fields can be bent or squeezed to an arbitrary extent 

^>j ' and without reflection using only dielectric structures. 

'NT ' 

o 

00 

o 
> 

X 



INTRODUCTION 



The coordinate transformation method (CTM) p| employs the invariance of Maxwell 
equations under coordinate transformations m to establish an equivalence between metric 

n n 

transformations and changes of material parameters [3| , [^ . The materials with parameters 
chosen to mimic a desired coordinate system are called the transformation media pj, |5|. 

CTM has been exploited for various computational problems such as: design of perfectly 
matched layers (PMLs) for simulations of open boundaries in finite-difference time-domain 

n Q n n 

and finite element methods 6[, l7j, [8j, [9|, simplifying the geometry of complex computa- 



tional domains 



io| 



M 



m 



representing waveguide bends and twists by equivalent 



straight segments IJ], [l5|, and formulating novel perturbation schemes for anisotropic 
materials [l6|. In l|], Pendry et al. suggested implementing transformation media as meta- 
materials and inspired the widespread use of CTM as an optical design tool. 

Several distinct strategies for the CTM-based design of metamaterial structures have been 
proposed. A continuous transforrnation of the whole space gives devices that are inhere ntly 
invisible (like invisibi 



Q 



ity cloaks [l|], electromagnetic field concentrators [17|], rotators 18 1 



[n 



and perfect lenses |19L yj). Various authors have previously applied CTM for design and 

Q, m m 



modeling waveguides [ij], [15[, [20j, |2l| . 
transformation method has been described 



22|. More recently, the embedded coordinate 
2311. where transformation media are embedded 
into surrounding space yielding devices that transfer the transformed fields from the devices 
to their exterior. 

In this paper, we apply CTM to domains bounded by transformation-invariant boundary 
conditions (BCs). Such an approach is inspired by the fact that in practice electromagnetic 
fields are always confined in space, like in waveguides and electromagnetic cavities or like 
electromagnetic beams. A refiectionless CTM-based device can have a twofold function: it 
can rearrange at will the field distribution within itself and it can yield an orthogonally 
transformed (rotated) outgoing wave. The beam squeezer is considered as a typical example 
of the former and the waveguide bend as the typical example of the latter. 

The CTM-based design offers the possibility of realizing a device with a given function 
(e.g. cloaking) in countless ways but it is a rule of thumb that the resulting prescription 
for material parameters is very hard to implement. For optical applications, probably the 
biggest problem is to obtain a material with a controlled permeability. Non-magnetic CTM 



devices can be designed, but only for TM waves. Therefore, as our main goal for the con- 
sidered structures (bend and squeezer) we choose achieving a non-magnetic realization. In 
Section III we show that using a reduced set of material parameters instead of the originally 
magnetic transformation media may be promising, provided the impedance-matching con- 
dition is reasonably satisfied. In Section IV we propose a different approach by considering 
transformations that yield non-magnetic transformation media. By giving examples of a 
non-magnetic waveguide bend and beam squeezer, we anticipate that in principle a struc- 
ture with arbitrary shape can be designed as non-magnetic, provided that its optical length 
is unrestrained. 

The numerical simulations in the paper have been done using the COMSOL Multiphysics 
FEM solver. 

THE COORDINATE TRANSFORMATION METHOD AND CONFINED DO- 
MAINS 

Figure [1] depicts a domain in (a) Cartesian coordinates x* {D') and (b) curvilinear coor- 
dinates x' (D) given by 

Domain boundaries are denoted by dD' and dD with BCs corresponding to either zero- 
field or perfect electric conductor (PEC), both of which are invariant under coordinate 
transformations. The electric permittivity and magnetic permeability of the medium within 
fl'a.e labeled w,th// and //,.e=pect,vely. 

The Maxwell equations are invariant under coordinate transformations [2|, whereas the 



material parameters in D read 



\iM 



e'^ = A''AlA],e' ^ fi'' = A-^^A^,^' ^ , (2) 

where A is the determinant of A*, and summation over repeated indices (running from 1 to 
3) is assumed. The electric and magnetic fields (one- forms) in D' are Ei' and Hi/, so the 
fields in D read 

E, = A^E,, H, = A^H,. (3) 
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FIG. 1: (Color online) (a) The domain D' in Cartesian coordinates x* , and (b) its image, the 
domain D in curvilinear coordinates x*. Blue lines stand for dD^ and dD. Transformation media 
given by Eq. ([2|) mimic the curvilinear coordinates x*. Red arrows show the propagation directions 
of an electromagnetic wave, assuming homogeneous e* ^ and /x* ^ . 

The above coordinate transformation can be physically implemented by interpreting x* as 
Cartesian coordinates and choosing the permittivity and permeability of the physical media 
within D (the transformation media) to coincide with e*-' and yU*-' ^. 

Figure [2] explains how transformed fields outside transformation media are obtained. The 
domain D' is divided into three subdomains, D[, D'^, D'^ and transformed to Di, D2, -D3, 
respectively. Within D[ and D'^ it is assumed that the transformation satisfies 



Ai A'j 6'' 



6'^\ 



A = 1, x'' eD[U D'^. 



(4) 



Equation (jl]) means that Di and ^3 have the same shape as D[ and D3, respectively. 
Subdomain D'2 is transformed to D2 so that it continuously connects Di and D3. Now Eqs. 
([2]) and (jlj) assert that £*•' and /i*-' for x* G -Di U -D3 are same as the corresponding e* ^ 
and /x*-' for x* G -D^ U 1^3 except that they may be rotated (together with D'^^ and/or D3). 
Since a continuous transformation ([T]) is established between entire D' and D, the fields 
are transformed according to ([3]) in Di and D3 in spite of the fact that media remained 
unchanged with respect to that in D'^ and D'.^. This way a recipe is established for a 
reflect ionless rotation of field distribution and polarization or change of propagation direction 
of a confined electromagnetic field that passes through a domain (-D2) with transformation 
media. 

The method of transferring modified fields away from transformation media described 



above is an a.. 
et al. I23J, 



ternative to the one based on embedded coordinate transformations by Rahm 



2J|. The benefit here is that the absence of reflection comes as a simple 



consequence of transformation continuity. 




FIG. 2: (Color online) Explaining how transformed fields are obtained outside transformation 
media. Dashed boundaries indicate that only relevant parts are shown while domains may extend 
arbitrarily beyond them. For simplicity, the media are assumed to be homogeneous in D' . (a) 
An electromagnetic wave (red arrow) propagates along a straight line in D' . (b) The wave is 
transformed by transformation media in D2 maintaining the orientation after passing through. 



Another application of the concept of bounded domains is in the design of structures with 
a desired external shape without affecting the way fields perceive the internal geometry. A 
typical example is in waveguide miniaturization 20|| or squeezing an electromagnetic beam, 
which we describe in detail below. 

The performance of any CTM-based device is, in principle, limited only by physical 
properties of the metamaterials implementing the transformation media. In the general case, 
parameters required for transformation media are such that causality implies the occurrence 
of both dispersion and absorption over any finite frequency band. These limitations have 
been reported for cloaking devices in [l|. Moreover, the fabrication is complicated since the 
metamaterial will need to be both anisotropic and with spatially varying properties in general 
case, as can be seen from ([21). Still, perhaps the biggest obstacle in optical applications is 
the difficulty in engineering the magnetic response. For these reasons, in CTM applications 
it is crucial to consider how can the material requirements be released or simplified even if 
it is sometimes achieved at the expense of departing from ideally required values. 

In the following two sections we consider two characteristic devices, the waveguide bend 
and beam squeezer, and discuss their implementation with respect to the above mentioned 
problems. Our main aim is to remove magnetism completely and that can be done for 
the TM polarized waves. In section III this is achieved by substituting the original (mag- 
netic) transformation media parameters with a reduced non-magnetic set having the same 
dispersion relation. In section IV we describe transformations that yield a non-magnetic 
transformation media. 



DEVICE DESIGN AND IMPLEMENTATION WITH A REDUCED SET OF MA- 
TERIAL PARAMETERS 

Waveguide bend 



The waveguide/beam bend shown in Fig. [3] has previously been treated in 2^, 2l|], 24| . 



251], 26|. The underlying transformation is 

X = y' sm{Kx') , y = y' cos {kx') 



(5) 



In this and all the following examples, we assume z = z' and free-space background e* ^ = 
eo^* ■' and /x* ^ = fioS^ ■' ■ By e^^ = fi^^ we will label the relative permittivity and permeability 
of the transformation media. 

D' C D = D' 
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FIG. 3: (Color online) (a) Domain D' from Fig. [2] (a) and (b) the corresponding domain D. The 
rectangular subdomain D2 is transformed to the annular segment shaped subdomain D2. The red 
arrows show the wave propagation direction when a metamaterial with parameters ([6]) is placed 
inside Z?2- 



The rectangle A'B'C'D' is transformed to the annular segment ABCD. L = \A'B'\ is the 
optical length of the bend while Ri = \0A\ and R2 = \0D\ are its inner and outer radius. 
If the bend angle is a, then kL = a. The relative permittivity is diagonal in cylindrical 
coordinates (r, 6, z) and reads 



e"^ = diag [{Kr) ^,Kr,{Kr) ^) , a, (3 = r,(j), z. 



(6) 



From now on we are concerned only with TM waves {z being the magnetic axis) for which 
only e"^^, e^'^ and fi^^ = e^^ are relevant. Following the procedure used in [271 and 28| for 
the cloak, a reduced parameter set for the bend from Fig. [3] (b) is found as 25|, 26|] 



[nr)' 



1, /x^^ = l. 



(7) 



In both the ideal ([6]) and reduced set of parameters ([7]), the free parameter k appears. 
Since parameters ([6]) are impedance-matched to vacuum irrespective of the value of k, it 
affects only the phase shift in passing through the bend. However, reduced parameters ([7]) 
are impedance- matched to vacuum only for /tr = 1, so the value of n determines the amount 
of reflection at the entrance AD and exit BC of the bend. The value of k that minimizes 
reflection depends on the incoming wave but for a general case, it is reasonable to set 

matching the impedance along the central line of the bend. Now we see that in case of 
reduced parameters, the bend geometry {Ri and R2) will play the deciding role in its per- 
formance through impedance mismatch away from r = Re- 

Figure H] shows numerical simulations results for a 90° waveguide bend. The second TM 
mode, TM2, with frequency 2GHz is excited at the left edge of structures. The simulated 
domains have PEC boundaries and are terminated by PMLs to simulate infinite waveguides. 
For reference. Fig. Il](al) and (a2) show the case where ideal parameters ([6]) have been used. 
The lower part of Fig. HI (bl) and (b2), shows the case of reduced parameters (JTj) together 
with ([H]). From the weak standing wave pattern in front of and the dominantly TM2 pattern 
behind the bend. Fig. H] (b2), we see that reflection and modal mixing are low, so the device 
works fine. 

A possible implementation of the non-magnetic bend with reduced set of parameters is 
as a structure comprised of many concentric annular layers with homogeneous and isotropic 
permittivities. The continuous variation of e^"^ given by (171) is first approximated by A^ 
annular layers with constant e"^"^ = {uri)^^ and thickness d, Ti being the inner radius of the 
layer and i = 1,2,. ..,N. To obtain the anisotropic permittivity, each of the A^ layers is 
further divided into n layers with homogeneous permittivities Ek and thickness 6kd, while 



X]fc=i ^fc = 1 i^k are the relative thicknesses), as in [29|. n is usually taken to be 2, but we 
consider the case of n = 3, as well. 

If the wavelength is large compared to d, each of A^ layers can be considered to be an 



anisotropic medium with effective dielectric permittivity 30|, 29| 



1 " X " 



k=i " k=i 



First we consider the case n = 2 and take that 61 = 62 = 0.5. From ([7]) we find Si and 62 




FIG. 4: (Color online) FEM simulation results for the 2GHz TM2 mode excited at the left edge 
of the structures: (a) the waveguide bend with ideal and (b) with reduced set of the material 
parameters. On the left side, (xl), the real part of magnetic field phasor is shown and on the 
right side, (x2), the magnetic field magnitude. The bend angle is 90°, Ri = O.lm, R2 = 0.3m and 
K = 1/0. 2m- in both cases. 



as 



£1 



[KTi 



62 = 1 + v/T^T^^. (10) 

If 1 — [nri)^"^ is to be positive within the bend (this corresponds to passive/lossless di- 
electrics), the condition ([8]) cannot be met. We have found that this is the main obstacle in 
realizing better bends with dielectric layers. 

Figure [5] (al) and (a2) show the simulation results for TM2 mode passing through a 
layered bend with n = 2 where k = 1/Ri {Ri = O.lm) and 40 layers (A^ = 20) have been 
used. Layer permittivities vary from 0.06 to 1.94. That the bending effect works very well 
is seen from a negligible modal mixing, as both the reflected and transmitted wave consist 
of almost only TM2. However, compared to the case shown in Figure |4] (bl) and (b2), there 
is a significant reflection (indicated by the pronounced standing wave in front of the bend) 
and the transmitted power is reduced. As we previously discussed, this is due to the choice 
of K. 

To improve on the implementation with ra = 2, we consider the case of ra = 3 anticipating 




FIG. 5: (Color online) FEM simulation results for the 2GHz TM2 mode excited at the left edge 
of structures, (a) The waveguide bend with k = 1/0. Im, A^ = 20 and 40 dielectric layers. Layer 
permittivities are gradually changed according to (jlOp . (b) The waveguide bend with k = 1/0. 08m, 
A^ = 20 and 60 dielectric layers. Layer permittivities are fixed (ei = 0.8, £2 = 0.05, and £3 = 8) 
while their thicknesses vary with r. (xl) show the real part of the magnetic field phasor while (x2) 
are the magnetic field magnitudes. 

that it will allow us to achieve spatially varying effective parameters ([7]) by varying 5k while 
keeping Sk fixed. The realization with fixed permittivities is clearly more favorable from a 
practical standpoint since it is easier to control the thickness of layers than their permittivity. 
Taking Sk as fixed parameters, from ([9]) and Ylik (^fe = 1, we find 5k as 



5i 

52 

53 

The solution for 5k for a given value of nr has to satisfy 

5k >0, k = 1,2,3, 
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(11) 



(12) 



because 5k represent the relative thicknesses so they must be positive. Equation (fT^ deter- 
mines which values of k can be obtained. An additional problem arises from the fact that 
Sk should also be positive. We have found that in case of negative Ek the effective medium 



9 



description ([9]) is ruined for propagation parallel to the layers both due to the extinction 
of waves in layers with negative permittivity and due to the excitation of surface states on 
their interfaces. 

It is straightforward to show that fITT]) . flT^ and e^ > can be simultaneously met only 
if Kr > 1. With (11 01) we have seen that (jH]) cannot be met in the case of n = 2, so now we 
see that it is not possible for n = 3 either. The same conclusion applies for any n. 

Simulation results for the bend with A^ = 20 and 60 layers are shown in Fig. [5] (bl) and 
(b2). The value of n is chosen to be 1/0. 08m. Layer permittivities used in simulation are 
£1 = 0.8, 62 = 0.05 and £3 = 8. The reflection is still higher than in Fig. H] (bl) and (b2), 
but lower than in the case of the n = 2 bend. The improvement over the n = 2 is probably 
due to using more layers. The main advantage of the n = 3 bend is, as noted before, in it 
being implemented using only three different homogeneous and isotropic materials. 

We have, so far, shown that this bend (the one discussed in literature so far) can, for 
TM polarization, be implemented using non-magnetic materials so that modal mixing is 
negligible, but with evident reflection depending on details of the TM wave. In Section IV 
we will further improve on this by considering a non-magnetic transformation media (for a 
different bend) that will remove reflection completely. 

Beam squeezer 

Here we use ideas from Section II to design a device for manipulating free-space elec- 
tromagnetic beams. Figure [6] shows the underlying transformation for a beam squeezing 
structure. The beam is compressed within the left trapezoidal subdomain, then it passes 
through the long channel, and flnally, it is expanded to the initial width within the right 
trapezoidal subdomain. Vacuum is assumed as the exterior of the device. Similar devices 
have independently been considered in [20|] and [Ml] except that stretching of the channel 
has not been investigated there. 

While the bend is a typical example of a CTM-based structure used for obtaining a 
rotated outgoing wave, here the main aim is to make a structure that complies with external 
requirements (width and length) leaving the outgoing wave unaffected. According to Fig. 
El two rectangles, the left (L), A'B'F'E', and the right one (R), C'D'H'G', are gradually 
narrowed down from \AE\ = \A'E'\ = wo to \BF\ = w and connected by the rectangle 
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FIG. 6: (Color online) Transformation underlying the beam squeezer, (a) Rectangular domain 
in Cartesian coordinates (a), and (b) the squeezed domain in curvilinear coordinates x*. The 
transformation media for the trapezoidal subdomains are ()14p and ()16p for the middle subdomain. 
The red arrows illustrate the beam propagation direction corresponding to the simulation in Fig. 

m 

between them (the channel) which has been narrowed down to w and stretched out from 
\B'C'\ = Iq to \BC\ = I. The side rectangles are transformed using 



X = x' + Xa, y-Yo = Aaiy' - Yq), A„ = a^x + fe„, a = R, L. 



(13) 



The introduced parameters are defined with a^ = {w—Wo)wq ^s' ^, &l = 1— a^x^, a^ = —a^, 
bR = 1- aRXD, XL = 0, xr = I - lo, s' = \A'B'\ = \C'D'\ = xb - xa = xd - xc, y = Yq is 
the symmetry axis of trapezoidal subdomains (in this case Yq = 0) while xp represents the 
X coordinate of the point denoted by P. The material parameters for the left and the right 
trapezoids are found as 






A-1 a^A-'iy - Yo) 

a^A-\y - Yo) alA'J{y - Y^f + A, 
A-i 



(14) 



The transformation of the rectangle B'C'G'F' to the rectangle BCGF is a simple rescaling 



x={x' -xb')^ + xb, iy-Yo) =-i{y' -Yq), 



(15) 



where ^ = 11^^ and 7 = wWq ^ are the stretching and compression parameters. The material 
parameters for the BCGF rectangle are 



e'^=li'^=dmg{ir\d,d 



(16) 
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The effect of geometrical transformations on electromagnetic fields can be discussed on 
the example of the channel in Fig. [71 Using flTSj) and ([3]) the field intensities in the channel 
are found to be 

Ex = ^^ Ex', Ey = 7~ Ey>, H^ = H^i, (17) 

with Hz unaltered since the z-axis is left invariant. This shows that the compression (stretch- 
ing) along the ^/-direction (x-direction) is accompanied by increased (decreased) field compo- 
nent intensity along that direction. The squeezing structure focuses the fields concentrating 
all the power flow to an arbitrarily small (not limited by diffraction) cross section. It could 
be used for various sensor applications where a strong signal improves the sensitivity or for 
imaging purposes. 

Figure [7] (al) and (a2) show the simulation of a Gaussian TM beam (source denoted by 
S) with oblique incidence passing through the squeezer with ideal material parameters. To 
better illustrate how the fields propagate through the structure, the compression and the 
stretching were set to moderate values, 7 = 0.15 and ^ = 39. The structure entrance is 
twice as wide as the beam, ensuring that the beam does not reach structure boundaries. 
In this way, the zero field BC is satisfied. Having passed through the narrow channel, the 
beam leaves the squeezer translated (due to stretching of the channel) along the x-axis 
retaining its initial propagation direction and the field distribution. Since the structure is 
theoretically perfect the slight reflection is only an artefact of numerical simulation and can 
be decreased with higher mesh density. However, it also indicates an increasing sensitivity to 
material parameter deviations (here occurring due to numerical discretization) with higher 
compression (small 7) or stretching (bigger ^). 

The reduced set of material parameters for the squeezer is obtained by setting fi^^^ = 1 
while changing the relevant components of permittivity so that the dispersion relation re- 



mains unaffected |27l |: 



XX _ zz XX xy _ zz xy yy _ zz yy /-|o\ 

From our experience with the bend, we expect that reflection may be significantly decreased 
by impedance-matching to vacuum along the central line of the device {y = Yq). As in the 
case of the bend, we find the impedance-matching condition to be that the optical length 
equals the physical length, i.e. ^ = I- 
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FIG. 7: (Color online) FEM simulation of the beam squeezer: (a) ideal parameters (7 = 0.15 and 
^ = 39) and (b) reduced parameters (7 = 0.15 and .^ = 1). (xl) Real part of the magnetic field 
phasor, (x2) magnetic field magnitude distribution. 5 is the source of the TM beam. The color 
maps are saturated for better contrast. 

Simulation results for the squeezer with reduced material parameters and ^ = 1 are given 
in Fig. [7] (bl) and (b2). As can be seen, the device retains its function but with a visible 
distortion of the transmitted beam. Also, appreciable reflection appears from the entrance 
and exit of the device, as indicated by standing wave patterns in front and within the device 
in Fig. [7] (b2). To highlight the overall field distribution, color maps have been saturated 
(regions with pronounced standing wave peaks have been cut off). In simulated cases with 
^ 7^ 1 (not shown) we have found strong reflection and obscured field pattern. Thus the 
proposed squeezer can be reahzed with reduced parameters, but it is crucial to have its 
impedance along the central line matched to the exterior. 

DEVICES BASED ON NON-MAGNETIC TRANSFORMATION MEDIA 

In the previous section we have attempted to substitute the ideal (magnetic) transforma- 
tion media parameters with a non-magnetic reduced set having the same dispersion a;(k). 
That approach yielded good results but there was always the problem of reflection due to 
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impedance mismatch. 

As in the previous section a free-space background is assumed, Sq^e^'^' = jIq^ji^'-'' = 5^'^' . 
As long as the transformation leaves the 2;-direction invariant, z = z', one of the optical axes 
of the transformation media is parallel to it and we have 

fi'' = A-\ (19) 

which is the only relevant permeability component for TM waves. Thus if 

A = 1, (20) 

a non-magnetic transformation media for TM waves is obtained. If we agree to evaluate 
the volume in x* coordinates as if they were Cartesian (see note after ([3]) in Section II) , the 
condition (I2U]) means that every subdomain of D' is mapped to a subdomain of D with the 
same volume. 

So far, this idea received little attention in the literature. Chen and Chan consider layered 



dielectric structures that act as non-magnetic transformation media in [32(]. In [33|, Luo et 
al. describe a cylindrical cloak with spatially invariant permeability /i^^ which is based on a 
similar idea, except that A = const > 1 because the cloaking transformation decreases the 
volume of a domain by creating a hole in it. 

We now show that the waveguide bend and beam squeezer can be designed with a transfor- 
mation that satisfies (!20l) and thus conceptually solve the problem of bending and squeezing 
TM beams using non-magnetic structures. 

Non-magnetic waveguide bend 



The TM beam shifter described in [23| is non-magnetic and (!20|) is satisfied. The shifter 
does not change a field propagation direction, it only translates the field. A change of a field 
propagation direction can be achieved by a shift of the polar angle in cylindrical coordinates. 
This is shown in Fig. [H] where the annular ring segment D'^ in Fig. [S] (a) is transformed to 
the subdomain D2 of the bent waveguide in Fig. [8] (b). The transformation reads 

where ^0 is the bend angle, /(r') an arbitrary continuous function of r', Ri and R2 are the 
inner and outer radius of D'2. The relative permittivity tensor e'i in cylindrical coordinates 

14 




FIG. 8: (Color online) (a) Straight waveguide structure is transformed to (b) a bent waveguide. 
Transformation of D'2 to D2 is given by (j2ip . (c) Curved layers within D2 for implementation using 
dielectric layers. 



[a, P = r,(f), z) is given by 



.a/3 



1 mr 

mr 1 + mr'^ 

1 



ec^ 



m 



df{r') 



f{R2)-f{Ri) dr> 



(22) 



for the field rotator. The bend 



The transformation (12T|) is the same as the one used in 
is obtained by cutting out a segment from the field rotator. This device can be implemented 
using curved isotropic and homogeneous dielectric layers when /(r') = In(r') 32|. It is 
shown in Fig. [H] (c) where alternating dielectric layers with permittivities e\ = 0.037 and 
£2 = 27.33 are assumed. 

Figure [9] (a) shows the simulation results for the non-magnetic waveguide bend with PEC 
boundaries. The 1.66GHz TM2 wave is excited on the left edge. The structure is theoretically 
ideal so the weak refiection is a numerical error that can be decreased by increasing the mesh 
density or by increasing the ratio R2/R1 so that the bending is less abrupt. Figure M (b) 
shows the simulation for a bend realized using 13 dielectric layers. There is a slightly higher 
refiection than in Fig. [9] (a). The bend in Fig. [9] (c) is realized with 26 layers giving 
practically the same results as the ideal structure from Fig. [9] (a). 

These results demonstrate that a refiectionless and non-magnetic waveguide bend for TM 
waves can be fabricated using only isotropic and homogeneous dielectric layers. Compared 
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FIG. 9: (Color online) FEM simulation results for non-magnetic waveguide bend: (a) ideal bend, 
(b) bend with 13 dielectric layers and (c) bend with 26 dielectric layers. Layer permittivities are 
constant, ei = 0.037 and £2 = 27.33. (xl) show the real part of the magnetic field phasor, while 
(x2) are the field magnitude distributions. 

to the bend with reduced parameters analyzed in Section III, the non-magnetic bend clearly 
shows a superior performance that comes at the expense of a more complicated geometry. 

Non-magnetic beam squeezer 

Figure [TU] (a) and (b) show the transformation used for the non-magnetic beam squeezer. 
The shaded rectangles in Fig. [10] (a) and corresponding shaded subdomains in Fig. [in] (b) 
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FIG. 10: (Color online) (a) Rectangular domain in x* transformed to (b) the domain in x* by 
transformation given in (|23|) and (|26|) . Few subdomains in (a) and (b) are shaded to emphasize 
that their areas are equal since A = 1. Simulation of a 7.5GHz TM Gaussian beam excited at S: 
(cl) Real part of the magnetic field phasor distribution, (c2) magnetic field amplitude distribution. 
Simulation parameters: 7 = 0.125, s' = 0.06m, \BF\ = 0.1m and / = 0.32m. 

having the same area, illustrate the point that A = 1. For rectangles A'B'F'E' (L) and 
C'D'H'G' (R) the transformation is given by 



X = — In 

(Ify 



x' + 



+ Ca, y = {a^x + h)y' , a = R,L, 



(23) 



with parameters aji = —ql = (1 — 7)/s', s' = \A'B'\, 7 = \BF\/\B'F', b^ = 1 — aLx'^^,, 
6ij = 1 — a^x'j^i, I = \BC\ and 



Cl = x^, In 

ttL V 1 - 7 



I , / s' \ , 2s' . / 7s' 



, Cr = Cl + 1 



In 



(24) 



(25) 



1 ~ 7 \1 "~ 7. 
The relative permittivity tensor of transformation media within ABEF and CDHG reads 

/-2 a^yf-' 

aayf-' aWf-^ + p 

1 

with fa = \aa\ exp {aa{x — Cq,)). 

The transformation of the rectangle B'C'G'F' to the rectangle BCGF is determined by 
the condition that the compression 7 is reciprocal to the stretching ^ = \BC\/\B'G'\: 

x = 7 X H , y = 'yy, (26) 

ol 7 
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so the parameters of media in the rectangle BCGF are 

£^^■ = diag(7-^7^1)• (27) 

The simulation of a TM Gaussian beam passing through the squeezer in Figure [10] (cl) 
and (c2) confirms the perfect performance of the proposed device. 

This completes the demonstration of non-magnetic realization of the two principal ex- 
amples in manipulating confined TM fields. The approach we have taken in this section is 
clearly better than the one with reduced material parameters that is widely used for the 
design of non-magnetic CTM-based devices. It might not be apparent immediately, but the 
condition A = 1 is quite lax and is equivalent to fixing the phase of the transformed field, 
which we had to do even in the case of reduced material parameters to minimize reflection 
(recall fixing k for the bend and ^ for the squeezer). 

In summary, the application of CTM to spatially confined electromagnetic fields was con- 
sidered. It was used to give various solutions for the waveguide bend and beam squeezer. 
Both approximate and exact non-magnetic realizations were found. In the case of structures 
with a reduced set of material parameters, the influence of impedance mismatch and its 
dependence on free parameters were investigated. Transformations with unit Jacobian de- 
terminant were recognized as promising for the design of non-magnetic devices and used to 
find novel non-magnetic solutions for the bend and squeezer. Several implementations using 
isotropic and homogeneous layered systems were considered. On the example of the bend 
with reduced parameters, it was shown that the inhomogeneity of effective permittivity ten- 
sor can be successfully controlled by varying only the layer thicknesses. The implementation 
of the non-magnetic bend with such layers was shown to be practically perfect. Its applica- 
tion is in optical waveguides that allow TM polarized waves. The non-magnetic squeezer is 
interesting for focusing and imaging applications since it is not limited by diffraction. 
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